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abstract 


This final technical report contains^hree 

deals with the p~^ t S 2 e dea i s e w ith the 3-D degenerated element, 
applications. Part 2 del the two major tasks that were 
These two parts constitute tne r i late d topic that was 

completed under the grant. *-4 rt a+-irtn is the development of 

initiated during the present investigation is the in 

a nonlinear material model. This topic is “11*1^^ and 

reference^are^includec^in each'part^^In^the^interest of^ 

Th^detail^an^additional'topics are described in the references 


cited. 
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PART 1 


GEOMETRICALLY NONLINEAR ANALYSIS OF LAMINATED SHELLS 
INCLUDING TRANSVERSE SHEAR STRAINS 

J. N. Reddy and K. Chandrashekhara* 

(A <LondtYU>zd veAiion o £ tkii papeA -li to app&OA in AIA A JouAnaZ , 19S4) 


SUMMARY 

The paper contains a description of a doubly curved shell finite element 
for geometrically nonlinear (in the von Karman sense) analysis of laminated 
(doubly-curved) composite shells. The element is based on an extension of the 
Sanders shell theory and accounts for the von Karman strains and transverse 
shear strains. The numerical accuracy and convergence characteristics of the 
element are further evaluated by comparing the present results for the bending 
of isotropic and orthotropic plates and shells with those available in the 
literature. The many numerical results presented here for the geomerti cal ly 
nonlinear analysis of laminated composite shells should serve as reference for 
future investigations. 

INTRODUCTION 

Laminated shells are finding increased application in aerospace, automo- 
bile and petrochemical industries. This is primarily due to the high stiff- 
ness to weight ratio, high strength to weight ratio, and less machining and 
maintenance costs associated with composite structures. However, the analysis 
of composite structures is more complicated when compared to metallic struc- 
tures, because laminated composite structures are anisotropic and character- 
ized by bending-stretching coupling. Further, the classical shell theories, 
which are based on the Kirchhoff-Love kinematic hypothesis (see Naghdi [1] and 
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Bert are known to yield deflections and stresses in Urinated 

are as L as 30, in error. This error is doe to tne neglect of transverse 

shear strains in the classical shell tneories. 

Refinements of the classical shell theories (e.9*. Love s fir ” 

mation theory [3], for shells to include transverse shear deformat, on 

Reissner [4-6]. Sanders [7] presented modified firs 
been presented by Reissner l j 

- <.< that removed an inconsistency (nonvanishing of 

second-approximation theories that removed 

small rigid-body rotations of the shell, existed in Love's first-approximation 

theory ^ ^ ^ ^ ^ „ la minated orthotropic composite shells is 
doe to Ambartsumyan [3,9]. m these works Ambartsumyan assumed that the in i 
vidual ortnotropic layers were oriented such that the principal axes^mat, 
rial symmetry coincided with the principal coordinates of the she 
surface. Dong, Pister, and Taylor [10] presented an extension of Oonn 
shallow shell theory [11] to thin laminated shells. Using the 
tegration of the elasticity equations, Widera and Chung [12 derive 
approximation theory for the unsy-etnc deformation of nonhomogeneous. an - 
tropic, cylindrical shells. This theory, when specialized to isotropic mate- 

rials reduces to Donnell's shell theory. 

joe effects of transverse shear deformation and thermal expansion roug 

the Shell thickness were considered by Zukas and Vinson [13], Dong and To 
[14] constructed a laminated orthotropic she,, theory that includes transverse 

This theory can be regarded as an extension of Love s 
shear deformation. This theory can u C 

first-approximation theory [3] for homogeneous isotropic shells. Other re- 
fined theories, specialized to anisotropic cylindrical shells, were presented 
p, Whitney and Sun [15], and Widera and Logan [16,17]. 
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The finite-element analysis of layered anisotropic shells, all of which 
are concerned with bending, stability, or vibration of shells, can be found in 
the works of Schmit and Monforton [18], Panda and Natarajan [19], Shivakumar 
and Krishna Murty [20], Rao [21], Siede and Chang [22], Hsu, Reddy, and Bert 
[23], Reddy [24], and Venkatesh and Rao [25,26]. Recently, Reddy [27] extend- 
ed the Sanders theory to account for the transverse shear strains, and pre- 
sented exact solutions for simply supported cross-ply laminated shells. All 
of these studies are limited small displacement theories and static analyses. 

In the present paper, an extension of the Sanders shell theory that ac- 
counts for the shear deformation and the von Karman strains in laminated an- 
isotropic shells is used to develop a displacement finite element model for 
the bending analysis of laminated composite shells. The accuracy of the ele- 
ment is evaluated by comparing the results obtained in the present study for 
isotropic and orthotropic plate and shell problems with those available in the 
literature. Numerical results for bending analysis of cylindrical and doubly- 
curved shells are presented, showing the effect of radius-to-thi ckness ratio, 

loading, and boundary conditions on the deflections and stresses. 

* 

A REVIEW OF THE GOVERNING EQUATIONS 
Consider a laminated shell constructed of a finite number of uniform- 
thickness orthotropic layers, oriented arbitrarily with respect to the shell 
coordinates (Si» 5 2 »C)» The orthogonal curvilinear coordinate system 
Ul»£ 2 »C) 1S chosen such that the and curves are lines of curvature on 
the midsurface c = 0. and e-curves are straight lines perpendicular to the sur- 
face c=0 (see Fig. 1). A line element of the shell is given by (see Reddy 
[27] 

(dS) 2 = [ r 1 + C/R^a^J 2 + [(1 + C/R 2 )a 2 d5 2 ] 2 + (dc ) 2 (1) 
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where a i and (i = 1 , 2 ) are the surface metrics and radii of curvature, 
respectively. In general, and are functions of \ ^ only. For the doubly 
curved shells considered in the present theory, and R^ are constant. 

The strain-displacement equations of the shear deformable theory of 
doubly-curved shells are given by 


e? + C«fi 


~ e 2 


e 4 = e 4 


o 
4 
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e 5 = e 5 


c 6 ’ c 6 * 


( 2 ) 
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(3) 


Here u^ denotes the displacements of the reference surface along £_.(5 = £) 

axes, and 4 ^ are the rotations of the transverse normals to the reference 
surface. In Love's first-approximation theories the parameter c Q is taken to 
be zero, and it is introduced only in the Sanders theory. 
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The stress-strain relat 

the form 


ions, transformed to the shelf coordinates, are of 

{») • CQMc) (4) 


where d' k) are the material properties of k-tn layer. 

The principle of virtual wort for the present problem is given by 

0 . \ ^ j (< 4 % * 4 % ♦ 4 % * * 4 k ’ 6E 5 

k>l C k _! 0 


- qSu-^a^d^d^dC 

= I [N.6c° t N 2 6e” ♦ N 6 c° ♦ * V‘2 + V'6 * Q l 6e 5 

Q 

+ Q 2 6eJ - q6u 3 ]a 1 a 2 d^ 1 ^ 2 


(5) 


(6) 


,„.u N anf i m. are the stress and 
where q is the distributed transverse load, ^ i 

moment resultants, and Q, is the shear force resultant: 


L 


(V N i) ■ J, { 

* 1 C K-1 
L 

Qi - I 

1 k-l C k-1 


c.[l,C)<lC . 1 = i ’ 2 * 6 


o.dc , i = 4,5, 


(7) 


K? = <2 = 

the 


where K, (i - 1.2) are the shear correction factors (taken to be S 
5/6), and are the ^coordinates of the k-th layer, and L :s th 

total number of layers in the laminated shell. 

u , s informative to note that the equations of equilibrium can be 

oerived from E„ (6, by integrate the displacement gradients in by parts 
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and setting the coefficients of 6u i (i = 1,2,3) and 6<i>. ('i = 1 , 2 ) to zero 
separately. We obtain [with c Q = ^ — -^-) and dx i = ct-jd^. ] 


aN l a Qi 

+ - c oV ^ 0 

* 3Np Qp 

3^ (N 6 + C 0 M 6 ) + 3^ + = ° 

8Q, SQn N. 


3M. dM, 

— - + — - -0=0 
dXj dx 2 y l 


dM, dM ? 

wf * wf • «2 ■ ° 


( 8 ) 


where 


«(“ 3 ) ■ 


5u 


ax. 


au. 


(N. — - + N, — i) + 

1 l ax. 6 ax' ax. 


au, 


au. 


( N 6ST + N 


2 ax. 


(9) 


The resultants (N^ , , Q i ) are related to (e°, <■) { i , j = 1,2,6) by 


N. = A. . e + B. 

1 ij J 1 J J 

"> * B u E ! + Vl 

Q 2 = A 44 e 4 + A 45 e 5 
Q 1 = A 45 e 4 + A 55 e 5 


(10) 

(ID 


Here A i j » B ij and (i.J = 1>2,6) denote the extensional , flexural- 
extensional coupling, and flexural stiffnesses of the laminate: 


C, 




k-l 


Q • j ) (l,C,C 2 )dC ( i » j - 1 , 2 , 6 ) 


( 12 ) 
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The boundary conditions. -rlv* «• *’«“•’ "* 

specifying either the essential boundary conditions (EBC) or the natura 

boundary conditions (NBC): 

EBC 

or 


Ur 


or 


NBC 

Vl + ( N 6 * C o M 6 )n 2 
N 2°2 + (N 6 + C o M 6^2 


au. 


u- 


or 


( N i ^7 )n l * (n 2 i^ 1 " 2 

au- 5u 3 

+ Cn — -In.. + (\ 77“J n 2 
^ 6 ax 2 J l 6 2>x 1 ^ 


+ Q 2 n 2 + q i n l 


or 


or 


Vi + M 6 n 2 


M 2 n 2 + Vl 


(13) 


.here denote the direction cosines of the unit normal on the boundary 

Of the midsurface of the shell. 

The exact for™ of the spatial variation of the solution of c,s. (81-U3). 
for the smal 1 -displacement theory, can be obtained under the folio., n g condi- 

tions (see Reddy [27]) laminates- i e., laminates 

( i) symmetric or antisymmetric cross-ply laminates. 

with 


a 16 - a 26 - b 16 - b 26 - o 16 - d 26 - a 45 - o. 

(ii) p r eely supported boundary conditions: 

N 1 (0,x 2 ) = N 1 (a,x 2 ) = M 1 (0 ,x 2 ) - '\( a,X 2 ) 
u 3 (0,x 2 ) - u 3 (a,x 2 ) - u 2 (0,x 2 ) - u 2 (a,x 2 ) - 0 


( 14 ) 
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V x l’°) = N 2 ( X r b) * M 2 (x 1» 0) = Mg( x x ,b) = 0 

u 3 ( x i,°) = u 3 (x 1§ b) = u^Xj.O) = Ujfxj.b) = 0 

<t 2 ( 0 ,x 2 ) = 4» 2 (a,x 2 ) = ♦ 1 (x 1 ,0) * ♦jtXj.b) = 0 ( 15 ) 

where a and b are the dimensions of the shell middle surface along 

the Xj and x 2 axes, respecti vely . The time variation of the load 
does not influence the spatial form of the solution. 

Note that the exact solution can be obtained only for cross-ply laminated 
shells with simply supported boundary conditions. For general lamination 
schemes, exact solutions are not available to date. 

FINITE-ELEMENT MODEL 

A typical finite element is a doubly-curved shell element in the x^x 2 - 
surface. Over the typical shell element the displacements 

(Ui, u 2 , u 3,4> 1 ,<t> 2 ) are interpolated by expressions of the form, 

N i 

u i = .1 u i 4 'j (x i * x 2 5 ’ 1 = 1 ’ 2 ‘ 3 

J ~1 

N * 

V J 0^j(x 1 ,x 2 ) , i - 1,2 (16) 

«J ^ 

where are the interpolation functions, and u 3 and are the nodal values 
of u^ and 4., respectively. For a linear isoparametric element (N = 4 ) this 
interpolation results in a stiffness matrix of order 20 by 20 . For a nine- 
node quadratic element the element stiffness matrix is of order 45 by 45. 

Substitution of Eq. ( 21 ) into the virtual work principle, Eq. ( 9 ) yields 
an element equation of the form 
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0(a)] (a) = {F} 


(17) 


, , r i U U T [k] the element stiffness matrix, 

where { 4 } - U-jt. ( u 2> • W > 1*1> • t*2 H ’ LKJ „ nf 

In the interest of brevity, the coefficients of 
and (F) is the force vector. In the interesu 

the stiffness matrices are inched 1" Appendix I. 

The element equations (17, can be assembled, boundary conditions can be 

imposed, and the results equations can be solved at each load step, 
that the stiffness matrix W it a function of the unknown solution vector 
U h therefore, an iterative solution procedure is required for each load 
step, in the present study, we used the direct iteration technique, which ca 

be expressed as 


[K((4} r )]W rtl ■ (Ft 


( 18 ) 


w „ ere (of denotes the solution vector obtained in the r-th iteration (at any 
given load step). At the beginning of the first load step, we assume that 
mo = (0} and obtain the linear solution at the end of the first iteration. 
The solution obtained at the end of the r-th iteration is used to compute the 
stiffness matrix for the (fl)-th iteration. At the end of each iteration 
(for any load step), the solutions obtained in two consecutive iterations 
compared to see if they are close enough to terminate the iteration and to 
in0 ve on to the next load step. The following convergence criterion is use in 

the present study: *> 

[ J Ia'J-aF* 1 ! 2 ' I |aF|’] 1/2 < 0.01 (19) 

L 1 1 i=l 


where N is the total number of unknown penalized displacements in the finite 
element mesh. 
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To accelerate the convergence, a weighted average of, the solution from 


last two iterations are used to compute the stiffness matrix: 

WyU) 1 " 1 ♦ U - r)(A} r )](4} r+1 ■ (F) (20) 

where y is the acceleration parameter, 0 < y < 1. In the present study a val- 
ue of 0.25 - 0.35 was used. 


NUMERICAL RESULTS 

Here we present numerical results for some sample problems. To illus- 
trate the accuracy of the present element, first few examples are taken from 
the literature on isotropic and orthotropic shells. Then results (i.., de- 
flections and strsses) for several laminated shell problems are presented. 

The results for laminated shells should serve as references for future inves- 
tigations . 

All of the results reported here were obtained using the double-precision 
arithmatic on an IBM 3081 processor. Most of the sample problems were an- 
alyzed using a 2 x 2 uniform mesh of the nine-node (quadratic) isoparametric 
rectangular element. 

1 . Bending of a simply supported plate strip (or, equivalently, a beam) under 
uniformly distributed load. 

The problem is mathematically one-dimen'sional and an analytical solution 
of the problem, based on the classical theory, can be found in Timoshenko and 
Womowsky-Krieger [28]. The plate length along the y-coordinate is assumed to 
be large compared to the width, and it is simply supported on edges parallel 
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to the y-axis . The following simply boundary cpoditions are 

assumed: + 

„ = » 2 = 0 along edges x = ± 127mm ' zl < 

AU inplane displacement degrees of freedom are restrained, ft o 

' me half Dlate is used to analyze the prob- 
four-node rectangular elements in the half plate 

+ , in Fia 2 The present result is m 
, era . The data and results are presented in Fig. 

good agreement with the analytical solution. 

2. rl aimed square ""' 1 - r ' in1forro l0a ^ _ ... 

^f^-^^^uadrant of the plate is modelled 

„,tn the 2 x 2 mesn of nlne-node elements (4 x 4 mesh of linear elements give 
almost the same result,. Pertinent data and results are presented in ig 
for side to thickness ratios a/, - TO and 500. The result for a/h ■ 500 
agreement with the results of way [29], The difference is attributed to the 
,. et that the present model Includes the Inplane displacement degrees of , ree- 

dom and transverse shear deformation. 

f1 g u re 4 contains transverse deflection versus load for Camped ortho- 

tropic. cross-ply, and ang, e-ply plates. The lamina properties are 


E = 25 x 10 4 N/mm 2 , E 2 


= 2 x 10 4 N/mm 2 , G 12 = G 13 = 1Q4 N/ ^ mm 


G 23 = 0.4 x 10 4 N/mm 2 , v l2 = 0.25. 

For the same total thickness the clamped orthotropic square plate is stiffer 
than noth two-layer angle-ply and cross-ply plates. 


13 





Deflection, -w/h (-w/h) 


Figure 3. Bending of clamped isotropic square plate under 
uniform load. 
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Deflection, -w (in mm) 


Figure 4. Bending of clamped orthotropic and laminated 
square plates under uniform load. 
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_ chell under po i>nt load^ 

3. »■?!» supported , i r r t l ” ,erica1 She1l ~ E - 

" ra of , he shel , is shown in Fig. 5. A uniform mesh of 
The pertinent data of the sneu 

. r thp pf f ect of three types o 
by 2 quadratic elements is used in a quadrant. The 

.. rontpr deflection and center normal stress 
simply supported conditions on the center deflect 

is investigated: 


SS-1 : u = w= 4>i = 0 


at y = b; v = w * 


= 0 


at x = a 


. . _ 0 at y - b; u - v - w - <t> 2 = 0 at x - a 

SS-2 : u = v = w - ♦ - u az J 


( 22 ) 


SS-3: v = w = 0 


at y « b; u = w = 


= 0 


at x 


Table 1 contains the results for the three Ooundary conditions. It is clear 
1 the results that al, three houndary conditions give virtually the same 
reS u,ts for a/h . ». and differ significantly (especially SS-1 ~ 

both SS-2 and SS-3) for a/h . 16. Thus, the effect is more in thict shells 

. tin c i c evaluated at point x 
. lie Thp stress o shown in Fi9* ^ 
than in thin shells. The stress a x 

y = 1.691" in the top layer 

^ ,i -i nHri r al shell linripr point load. 

4. Simply supp orted isotropic cy Vj 

The geometry and finite-elementmesh of the she,, are shown in F,g 

Once again, the effect of various simply supported boundary conditions < 

a stresses for the problem is investigated using a uniform 
tne deflections and stresses 

Th rocnits are presented in Table 2. 
mpsh of 2 X 2 quadratic elements. The resul P 

a here <R - 2540 a . 254 . h - 12 .7 ) , the boundary 

the geometry and loading used here (R - 254U, 

hionificant effect on the solution. Boundary conditions 
conditions have very significant eriecc 

SS-2 and SS-3 give almost the same results whereas SS-1 gives about 2-1/2 
times the deflection given by SS-2 or SS-3 boundary conditions. 
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Table 1. Effect of various simply supported boundary conditions on the center 
deflections and normal stress in spherical shells under point load 
(E = 10 7 psi , v = 0.3) . 


Load 

P/h 2 

Solution 

SS- 

a/h=160 

-1 

a/h=16 

SS- 

a/h=160 

2 

a/h=16 

SS- 

a/h=160 

■3 

a/ h= 16 

4,000 

-w* 

0.0155 


0.0152 

- 

0.0152 

- 


-a * 

893 

- 

984 

- 

894 

- 

8,000 

X 

-w 

0.0329 

0.0349 

0.0324 

0.0255 

0.0324 

0.0258 


-a v 

1,880 

6,535 

1,882 

6,015 

1,882 

6,031 

12,000 

X 

-w 

0.0529 

- 

0.0522 

- 

0.0521 

- 


-a 

2,980 

- 

2,985 

- 

2,986 

- 

16 ,000 

X 

-w 

0.0760 

0.0793 

0.0752 

0.0520 

0.0751 

0.0525 


-a v 

4,220 

13,230 

4.228 

12,200 

4,229 

12,240 

20,000 

X 

-w 

0.1038 

- 

0.1028 

- 

0.1027 

- 


-a 

5,657 

- 

5,671 

- 

5,672 

- 

24,000 

X 

-w 

0.1364 

0.1083 

0.1354 

0.0792 

0.1353 

0.0800 



7,268 

20,110 

7,289 

18,500 

7,291 

18,550 

23,000 

X 

-w 

0.1761 

- 

0.1752 

- 

0.1751 

- 


-a v 

9,128 

- 

9,160 

- 

9,162 

- 

32,000 

-w x 

0.2234 

0.1472 

0.2227 

0.1072 

0.2227 

0.1083 


‘ a x 

11,180 

27,170 

11,220 

24,930 

11,230 

25,000 

* w( 0 ,0) 

. a x ( A, A); 

A = 1.691 







Table 2. 

Effect of various types of simply supported boundary conditions on 
the deflections and stresses of anisotropic cylindrical shell under 
point load. 

Load ,P 

SS4 

1 

SS-2 


SS- 

3 

(N) 

-w(mm) 

-o ( N/mm 2 ) 

-W 

’°y 

-w 

"°y 

250 

2.5804(2) 

2.868 

0.6544(4) 

1.706 

0.6698(4) 

1.706 

500 

5.1626(2) 

5.713 

1.3533(4) 

3.478 

1.3843(4) 

3.477 

750 

7.7343(2) 

8.506 

2.1057(4) 

5.327 

2.1522(4) 

5.321 

1,000 

10.278(2) 

11 .210 

2.9234(4) 

7.265 

2.9855(4) 

7.242 

1,250 

12.733(2) 

13.80 

3.8241(4) 

9.312 

3.9017(4) 

9.288 

1,500 

15.204(2) 

16.25 

4.8349(4) 

11.50 

4.9279(4) 

11.46 

1,750 

17.560(2) 

18.560 

6.0331(5) 

13.91 

6.1423(5) 

13.85 

2,000 

19.843(2) 

20 .730 

7.5316(6) r 

16 .66 

7.6610(6) 

16.57 
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Figure 5. 


nding of a simply supported (SS-3), 
otropic, spherical shell under point 

ad. 


(psi) 
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Figure 6. Geometry of a cylindrical shell. 
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6 -i, t .h HotroEic C jrUndrical shell under uniform loadi ga^. 

figure 7 certains the pertinent data and resuits for a clawed cylindri- 
cal shell (isotropic) subjected to uniform load. The results are compared 
eith those obtained by Dhatt [30]. The agreement is very good. 

6 . n amneri orthotropic cylindrical shell -lihjerrpd to internal pressu re. 

Figure 8 contains the geometry and plots of center deflection and center 
stress versus the interna, pressure for the problem. The orthotropic materia, 
properties used in the present study are: 

r » 7.5 x 10‘ psi. E- * 2 x 10‘ psi. 0j 2 ■ 6 U - G 23 - 1-25 * 10« psi 


V,„ = 0.25 


(23) 
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. thp o x 2 mesh of quadratic elements, is in 

The present result, obtained using the 2 x 2 mesn 4 

t- uitn that obtained by Chang and Sawamiphakdi [31]. 
excellent agreement with that ooumeu j 


is used in the analysis (with SS-3 boundary 


7 . .,.., r r X0V90 V 0! ~/gl cross^ -pP-riral shell subjected _to 

uniformly distributed load^ 

The following geometrical data 
conditions) : 

R . r 2 , r = 1,000 in., a = b = 10O in., h - 1 in. (24) 

Individual layers are assumed to be of egual thianess (h, - W»>. with the 
zero-degree layers being the inner and outer layers. The following two sets 
of orthotropic-materia, constants, typical high modulus graphite epoxy mater, 
a, (the ratios are more pertinent here), for individual layers are used: 
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3 4 5 6 -Oy(N/mm 2 ) 

Deflection and Stress 


Figure 7. Bending of a clamped, isotropic, cylindrical shell under 
uniform load. 
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Mat.-l: E l = 25 x 10 6 psi, E 2 = 10 6 psi , G 12 = G 13 » 0.5 x 10 6 psi 

G 23 = 0.2 x 10 6 psi, v 12 = 0.25 (25) 

Mat. -2: E 1 = 40 x 10 6 psi, E 2 = 10 6 psi, G 12 = G 13 * 0.6 x 10 6 psi 

G 23 = 0.5 x 10 6 psi, v 12 = 0.25 (26) 

Figure 9 contains plots of center deflecton (w/h) versus the load parameter 
(P = q 0 R 2 /E 2 h 2 ) for the two materials. Shell constructed of Material 1 
deflects more, for a given load, than the shell laminated of Material 2 
(because Material 2 is stiffer), and consequently experiences greater degree 
of nonl i neari ty . Note that the difference between the nonlinear deflections 
of the two shells increase nonlinearly, indicating that the shell made of 
Material 2 can take much more (ultimate) load than apparent from the ratio of 
moduli of the two materials, E^/E^. 

8 . Effect of various simply-supported boundary conditions on the deflections 
of two-layer cross-ply spherical shells under uniform load. 

As pointed out in Problems 3 and 4, the transverse deflection is sensi- 
tive to the boundary conditions on the inplane displacements of simply sup- 
ported shells. To further illustrate this effect for laminated shells, a set 
of four types of boundary conditions are used, and the results are presented 
in Table 3. Here SS-4 has the following meaning: 


SS-4 


w = = 0 on x = a 

m = $ = 0 on y = b 


(27) 
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Figure 9. Bending of nine-layer cross-ply 
[0V90V0 0 /. . .] spherical shell 
subjected to uniformly distributed 
load. 
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Table 3. Effect of various simply-supported boundary conditions on the trans- 
verse deflections of cross-ply [0°/90°] spherical shells under 
uniform load (Material 1; shell dimensions are the same as those in Fig.* 


a 


-w 1 

(in.) 


H o 

(psi ) 

SS-1 

SS-2 

SS-3 

SS-4 

0.50 

0.3344 

0.04246 

0.04257 

0.4592 

0.75 

0.5757 

0.06599 

0.06617 

0.8255 

1.00 

0.9485 

0.09144 

0.09171 

1.3845 

1.25 

1.6529 

0.11926 

0.11966 

1.9589 

1.50 

2.2826 

0.15008 

0.15063 

2.3597 

1.75 

2.6421 

0.18478 

0.18556 

2.5951 

2.00 

2 .8499 

0.22473 

0.22584 

2.8074 

2.25 

3.0764 

0.27425 

0.27593 

3.0284 

2.50 

3.2432 

0.33534 

0.33795 

3.1948 

2.75 

3.4214 

0.42970 

0.43487 

3.3719 


26 


. ^ ? and SS-3 give almost the same deflections. 

Once again we note that SS-Z ara a 

„ , and SS-4 give deflections an order of magnitude hig 
Boundary conditions SS-1 and 55 4 g 

K ? and SS-3 Thus, boundary conditions SS-2 and 55- 

er than those given by SS-2 and S5 5. '"«»• 

make the shell quite stiffer. 

_ . , r ac o /a^ 0 1 simp 1 y-supported_ 

9 . cross- ply fOVWI and angle-ply [-45 /*Lii P / 

( spherical shells^ 

Figure 10 contains the pertinent data and results (vntti different sea 
for the cross-ply and angle-ply snells (of Material 2). » ^ interesting to 
note that the type of nonlinearity exhibited by the two shells is quite di 
ferent; the cross-ply shell gets softer whereas the angle-ply shell gets 

stretching coupling due to the lamination scheme (B 22 - - * n nonzero for the 
cross-ply shell and B 1S and „ 26 are nonzero for the angle-ply shell), the 
angle-ply experiences shear coupling that stiffens the spherical shell re a- 
t1vel y more than the norma, coupling (note that, in genera,, shells get softer 

under externally applied inward load). 

Figure 11 contains plots of center deflection, normal stress and 

Shear stress (u y z ) at x . y - 5.283" versus load for two-layer cross-ply 
(0°/90°) spherical shell (Material 1) uhder poiht load at the center of the 
shell. The nonlinearity exhibited by the stresses (especially c yz ) 1ess 
compared to that exhibited by the transverse deflection. 


Figures It ana u ' y 

cross-ply [0°/90°] and angle-ply [-«•/«*] damped cylindrical shells under 

uniform load. The load-deflection curve for the cros-ply she,, resembles that 
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Deflection, -w (10'^in) * angle-ply 



0.0 0.1 0.2 0.3 0.4 0,5 cross-ply 


Deflection, w (in in.) 


Figure 10. Bending of two-layer cross-ply and angle-ply, 
simply supported (SS-3) spherical shells under 
uniform load. 




Figure 11. 


Bending of a cross-ply (0°/90 o l 
(SS-3, Material 1), under point 
for the shell dimensions) 


spherical shell 
load. (see Fig. 10 
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Figure 12. Bending of a clamped angle-ply [-45V45 0 ] cylindrical 
shell under uniform load. (Material 1) 
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Figure 13. Bending of a clamped cross-ply [0°/90°l cylindrical 
shell under uniform load (Material l] 



of the isotropic shell in Fig. 7, but exhibits greater de-gree of nonlinearity 
(being stiffer). The angle-ply shell exhibits different type of nonlinearity 
(softening type) for all loads. 


11 . Quasi-isotropic, clamped, cylindrical shell under uniform load. 

Two types of quasi-isotropic clamped cylindrical shells are analyzed: 


Type 1: [0°/45 o /90 o /-45 o ] 

Type 2: [0°/+45°/90] synu 


(28) 


Material 1 properties are assumed for each lamina (8 layers). The geometric 
data and results are presented in Fig. 14. Compared to the results presented 
in Figs. 12 and 13, the quasi-isotropic shells have the 'near-inflection' 
point at higher loads; the load-deflection curve has essentially the same form 
as that of the cross-ply shell (see Fig. 12). 


CONCLUSIONS 

A shear-flexible finite element based on the shear deformation version of 
the Sanders' theory and the von Karman strains is developed, and its applica- 
tion to isotropic, orthotropic, and laminated (cross-ply and angle-ply) shells 
is illustrated via numerous sample problems. Many of the results, especially 
those of laminated shells, are not available in the literature and therefore 
should serve as references for future investigations. From the numerical com- 
putations it is observed that boundary conditions on the inplane displacements 
have significant effect on the shell deflections and stresses. Also, it is 
noted that the form of nonlinearities exhibited by different lamination schemes. 
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Figure 14- 


Bending of clamped quasi -isotropic cylindrical shells 
nnHpr uniform load. 
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APPENDIX I 


St i ff ness Coefficients/. 


au- 


1 au 3 


, „ f 1 u 3 f , a _ 

Let f i = i d^rr * 2 2 ax 2 


CK U] = A U [SH] * A 16 (CS»] * [S 2 ‘]) * a 66 CS»] 

- C 0 (B 16 ([S 1J ] + CS 2 ‘]) * 2B 66 [S 22 3 - c o 0 66 [S“]) + — CS 0 °] 

[ K 12] - A U [S‘2] + A 16 [S“] * A 26 CS 22 ] + A 66 [S 2 ‘] 

- c o( B 26 [S22] • B 16 CSll] + c 0°66 [s21]) + Vl CS ° 0:i 


[K 13] , f 1 (A u CSll] ♦ A 16 [$ 12 ] + CS 21 ]) + A 66 [s22;l) 

♦ f 2 (A 12 [S 12 ] + A 16 [SU] ♦ A 26 [S 22 ] + a 66 [s 21 :) 

+ i- (A U [S 10 ] + a 16 [s 2 »]) + ^ (A 1Z CS‘»] + A 26 [S 20 ]) 

. c (!l« [S »] * £ CS 2 °]) - t- (A^ * *55^ 

o'-R^ ^2 1 

- c 0 [VB 16 [S 2 ‘] * B 66 [S 22 ]) ♦ f 2 (B 26 ts 22 : * B 66 [S 22 D] 

[Kl‘] . B n [S“] +b 16 ([s> 2 ] * cs 2 ']) * b 66 [s 22 3 

- C 0 (D 16 [S 2 >] + o 66 rs 22 ]) - A 55 cs«»] 
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[K‘S] ■ 6 12 [Sl2] + B 16 [S“] + B 26 [S 22 ] ♦ BggCS 21 ] 

" c of°26^ S22 ^ * °66^ S21 ^ ' R^ A 45^ S °°^ 

[K 2 ‘] = [K‘ 2 ] T 

[K 22 ] = A 22 [S 22 ] * A 26 ([S‘ 2 ] * [S 2 ‘]) * A 66 [Sl>] * 2c 0 B 66 CSl»] 

* c o (B 26 ([Sl 2 ] * CS 21 ]) ♦ C 0 D 66 [Sll]) -^-[s»»] 

R 2 

[K 22 ] = f !(A 12 CS21] + A 26 [S 22 ] ♦ A 16 [S“] ♦ AggtS 12 ]) 

* f 2 (A 22 [s 22 ] ♦ a 26 ([s 2 ‘] * [si 2 ]) + a 66 [s"]) 

* «: ( V s2 °3 * A 16 [Sl0] ) + R^ ( A 22 [s20;l + A 26 CSl0] > 

* ‘o (sf ♦ sf )ts‘ 0 ] Rj ( A 44^° 2 ^ + A 45 CS 01 ]} 

* c„[VB 16 [Sn] * B 66 [S> 2 ]) * f 2 (B 26 [Sl 2 : * B 66 [S“])] 
[K 24 ] - B 12 CS 21 ] * 3 26 [S 22 ] * B 16 [S“] + B 66 [S> 2 ] 

+ W 5 "] + “k» 1J 3) * ^ 

[K 25 ] = B 22 [$ 22 ] ♦ B 26 ([S 2 ‘] ♦ [$> 2 ]) + B 66 [S“] 

* * °66 [Sll] ) -^ A 4 4 C S0 °] 
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U 31 3 NL ■ C 2 • 


NL = Nonlinear portion 


of the matrix 


[k 32 0 nl - tZK 23 3 ^L 


[H 33 ] * A^CS 32 ] * A 55 CSU: + ^ * A45CS ] 

* ZtS'LHV? * 2A W f l f 2 + 

♦ 2([S 12 ] * [S 21 ])[f^ W * ( A 12 * A 66 )f l f 2 * 


+ 2[S 22 ](A 66 f5 * 2A 26 f l f 2 + A 22 f 2 5 

A,, A 


, A,, a 12 , .^2 + 222 -n 

* [S°°][^ (-^ * -R^l R 2 ( «! R 2 n 

k n A 12 f rh§- + ^ 2 §oi 

+ f[s°i] + 2 [s io 3)[v4: + “r; ) + 2( R i r 2 


Nfi A 26 . + 

+ (CS 02] + 2[S20])[f 1 (^- ♦ -q) + V Ri fy J 


[K 34 ] = a 55 [s 10 1 + A 45 cs 2 °] 

* 2 f 1 (B ll [s“] * B 16 ([S 12 3 * ts 21 3) * 8 6« cs22;n 

+ 2f 2 (B 12 CS 21 ] ♦ B 66 [S 12 3 + B 26 [S 22 3 + B 16 [S 
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DC”] = a 45 [s 10 ] + a 44 [s 20 ] 


- 2f 1 (B i2 [Si2] + B 66 [ S 21] + B 16 [S“] + B 26 [S22]) 


+ 2f 2 (B 22 CS22] + CS12] + B 9fi [S2l] + B^CSll]) 
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26 1 


66 


B 16 . 8 26 ■ 


B 12 a B 22 • 




[K 41 ] = [Kl 4 ] 1 , [K 42 ] = [ K 24 ] T , [K 4 3 ] nl = \ CK3 4 ]J l 


[K 44 ] « D n [S“] + D lfi ([Sl2] + [S21]) + D fifi [S22] + A„[SOO] 


; 16 C 


66 l 


55 L 


[K 45 ] = d 12 [s 12 3 + O^CS 11 ] + D 9A [S22] + D«[S21] + A^CS”] 


'16 


26 L 


66 


45 


[K 5 1] = [K 15 ] T , [K 52 ] = [<25]^ [K 53 ] nl = [K 35 ]J l 


[K 54 ] = [K 45 ] 


[KS5] = d 22 [S22] + 0 ?fi ([Sl2] + [S21]) + D fifi [S“] + A Ad [S00] 


26 


66 1 
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f K a Pl = rK Pa 1 T 

LKl J Linear LK J| < 


Linear 


where 


S U ' I e — ^ dX l dX 2 • S ij * / , *1*J d, l dx 2 

n e ^ P ' r\ ® 
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u should be noted that although fj and f 2 are shown factored outside the 
matrices, in the evaluation of the coefficients by the Gauss quadrature fj and 
f 2 are considered as parts of the integrals. For example f^CS"] is 

evaluated by ^ gy 

/ f 1 A u P i P j dx 1 dx 2 i | £ J i A U [(^)* i *j] Xl =Z I .x 2 .Z J “l' , J detJ o 
Q 6 

where N is the number of Gauss points, Wj and VI, are the Gauss weights, Z, ahd 

Z, are the Gauss points, and 0 o is the Jacobian of the transformation. 

J 


WPtJNRKCl 
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PART 2 


ANALYSIS OF LAMINATED COMPOSITE SHELLS 
USING A DEGENERATED 3-D ELEMENT 

W. C. Chao* and J. N. Reddy 
Department of Engineering Science and Mechanics 

[Ttvii papa -cl to appear In Int . Jou/maZ o£ Numaical Methods in Engng . ) 

SUMMARY 

A special three-dimensional element based on the total Lagranglan 
description of the motion of a layered anisotropic composite medium is 
developed, validated, and employed to analyze laminated anisotropic 
composite shells. The element contains the following features: 
geometric nonlinearity, dynamic (transient) behavior, and arbitrary 
lamination scheme and lamina properties. Numerical results of nonlinear 
bending, natural vibration, and transient response are presented to 
Illustrate the capabilities of the element. 

INTRODUCTION 

Composite materials and reinforced plastics are increasingly used 
In automobiles, space vehicles, and pressure vessels. With the increased 
use of fiber-reinforced composites as structural elements, studies 
involving the thermomechanical behavior of shell components made of 
composites are receiving considerable attention. Functional 
requirements and economic considerations of design have forced designers 
to use accurate but economical methods of determining stresses, natural 
frequencies, buckling loads etc. 


Graduate Research Assistant; presently at the University of Dayton 
Research Institute 
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Majority of the research papers In the open literature on shells 
is concerned with bending, vibration, and buckling of isotropic 
shells. As composites materials are making their way into many 
engineering structures, analyses of shells made of such materials 
becomes important. The application of advanced fiber composites in jet 
engine fan or compressor blades and high performance aircraft require 
studies involving transient response of composite shell structures to 
assess the capability of these materials under dynamic loads. 

Finite-element analysis of shell structures in the past have used 
one of the three types of elements: 1. a 2-D element based on a two- 

dimensional shell theory, 2. a 3-D element based on three-dimensional 
elasticity theory of shells; and 3. a 3-D degenerated element derived 
from the 3-D elasticity theory of shells. The 2-D shell theory is 
derived form the three dimensional continuum field equations via 
simplifying assumptions. The simplifications require the Introduction 
of the static and kinematic resultants, which are used to describe the 
equations of motion. The unavailability of a convenient general 
nonlinear 2-D shell theory makes the 2-D shell element restrictive in 
its use. The degree of geometric nonlinearity included in the 2-D shell 
element is that of the von Karman plate theory. In contrast to the 2-D 
shell theory, no specific shell theory is employed in the 3-D 
degenerated element; instead, the geometry and the displacement fields 
are directly discretized and interpolated as in the analysis of 
continuum problems. 

Finite-element analyses of the large-displacement theory of solids 
are based on the principle of virtual work or the associated principle 
of stationary potential energy. Horrigmoe and Bergan [11 presented 
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classical variational principles for nonlienar problems by considering 
incremental deformations of a continum. A survey of various principles 
in incremental form is presented by Wunderlich [2]. Stricklin et al . 

[3] presented a survey of various formulations and solution procedures 
for nonlinear static and dynamic structural analysis. The formulations 
Include the pseudo force method, the total Lagranglan method, the 
updated Lagrangian method, and the convected coordinate method. 

The only large-deflection analyses of laminated composite shells 
that can be found in the literature are the static analysis of Noor and 
Hartley [4] and Chang and Sawamiphakdl [5]. Noor and Hartley employed 
the shallow shell theory with transverse shear strains and geometric 
nonlinearities to develop triangular and quadrilateral finite 
elements. Chang and Sawamiphakdi presented a formulation of the 3-D 
degenerated element for geometrically nonlinear bending analysis of 
laminated composite shells. The formulation is based on the updated 
Lagranglan description and it does not include any numerical results for 
laminated shells. 

From the review of the literature It Is clear that first, there 
does not exist any finite-element analysis of geometrically nonlinear 
transient response of laminated anisotropic shells, and second, the 3-D 
degenerated element is not exploited for geometrically nonlinear 
analysis of laminated anisotropic shells. In view of these 
observations, the present study was undertaken to develop a finite- 
element analysis capability for the static and dynamic analysis of 
geometrically nonlinear theory of laminated anisotropic shells. A 3-D 
degenerated element with total Lagrangian description is developed and 
used to analyze various shell problems. 
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The primary objective of this section Is to review the formation 
of equations governing geometrically nonlinear motion of a continuous 
medium. In the Interest of brevity only necessary equations are 
presented. For additional details the reader Is referred to References 

( 6 - 101 . 

We describe the motion of a continuous body In a cartesian 
coordinate system. The simultaneous position of all material points 
(l.e., the configuration) of the body at time t Is denoted by C t , 
and 0 and denote the configurations at reference time t - 0 and 

time t ♦ at. respectively (see Fig. 1). In the updated Lagrangian 
description all kinetic and kinematic variables are referred to the _ 
current configuration at each time and load step. In the total 
Lagrangian description all dependent variables are referred to the 
reference configuration. The updated Lagrangian is more suitable for 
motions that involve very large distortions of the body (e.g.. high- 
velocity impact). The total Lagrangian Is more convenient for motions 
that involve only moderately large deformations. In the present study 

the total Lagrangian formulation is adopted. 

„ ere we present a derivation of the equilibrium equations at 

different time steps using the total Lagrangian approach.^ The ^ 

coordinates of a typical point In C t Is denoted by l x - ( - r *3>- 
The displacement of a particle at time t Is given by 


l u . l x - °x or t u ( « *x ( -°x, 


( 1 ) 


The increment of displacement during 

t+At 


time t to t + it Is defined by 


u i = 


u. - u. 


( 2 ) 
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The principle of virtual displacements can be 'employed to write the 
equilibrium equations at any fixed time t. The principle, applied to 
the large-displacements case, can be expressed mathematically as 


J 


t+At;; 


u, tu, dV 0 ♦ s 


«( t+4t Mj) dV o 


- s 


t+At 


T, dA 0 ♦ ; v 0 0 t+lt F, ««, dV 0 
0 


(3) 


where summation on repeated Indices Is implied; Y 0 . A Qt and p q denote, 
respectively, a volume element, area element, and density in the initial 
configuration, S^- are the components of second Piola-Kirchhoff stress 
tensor, z.. the components of Green-Lagrangian straintensor, T i the 
components of boundary stresses, and F i are the components of the body 
force vector; the superposed dots on Ui denotes differentiation with 
respect to time, and 6 denotes the variational symbol. In writing Eq. 
(3) it is assumed that e . . is related to the displacement components by 

the kinematic relations 


t+At 


1 /t+At . t+At + t+At 

“i.j* “J.i "• 1 


t+At 


Vj> 


(4) 


where u, • = 3U./ax.. The strain components can be expressed in 

1 * 3 i j 

terms of current strain and incremental strain components as 


t+At _ 1 /t.. . tii + ^u .) 

e ij ‘ 2 ( U i,j «.i 

+ \ (u i , j + u j,i + tu m,i U m,j + U m,1 Nn.j* + ? Vi Vj 




( 5 ) 
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where e^ and denote the linear and nonlinear incremental strains. 
The stress components t+At S^ can be decomposed into two parts: 


t+At, 


*1J 


S ij + S ij 


( 6 ) 


where S^- is the incremental stress tensor. The incremental stress 
components are related to the incremental Green-Lagrange strain 
components, = e^ + by the generalized Hooke's law: 


S 1j * C i jkz c ki* 

where C 1 are the components of the elasticity tensor. Using Eq. (4)- 
(7), Eq. (3) can be expressed in the alternate form 


L °0 t+lt V U 1 dv o + ! + 1 ’k,‘«ij>< iv » 

0 0 

+ J V **11 “u *0 ■ * - ^ ts lj*"u«o < 8 > 

o w o 

where 6W is the virtual work due to external loads. 

FINITE-ELEMENT MODEL 

Geometry of the Element 

Consider the solid three-dimensional element shown in Fig. 2. 

The coordinates of a typical point in the element can be written as 


x, = 


>1 




its. ( *4 


l-c 


2 <Vtop + j f 1 *j (5 l- 5 2> 


^ x i ^bottom 


( 9 ) 


where n is the number of nodes, are the finite-element 

interpolation (or shape) functions, which take in the element, the value 
of unity at node i and zero at all other nodes, ^ and are the 
normalized curvilinear coordinates in the middle plane of the shell, 
and ; is a linear coordinate in the thickness direction and xj, xl, and 
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are the global coordinates at node i. Here and t are assumed 

to vary between -1 and +1. Now let (see Fig. 2) 


'3k 


i< 


«*k»top - W) 


k 7 bottom 


?3 * V3/IV3I 


( 10 ) 


where vl^ is the k-th component of the vector v^. Then Eq. (9) becomes 

X 1 = 1 [ *j (x i ) m 1 d + "jl h J*31 l <"> 

J 

where hj is the thickness of the element at node j. For small 
deformation, the displacement of every point in the element can be 
written as 


*j |u i + ' 7 <*11 s 2 - 'si 9 ^ 1 < 12 > 

where e-j and e 2 are the rotations about (local) unit vectors e^ and e^, 
respectively, u-i, u 2 , and u 2 are the displacement components 
corresponding to the global coordinates x-j, x 2 , X 3 directions 
respectively, and u j , u 2 and u^ are the values of the displacements 
(referred to x) at node i. In writing Eq. (12), we assumed that a line 
that is straight and normal to the middle surface before deformation is 
still straight but not necessarily ‘normal 1 to the middle surface after 
deformation. The strain energy corresponding to stress perpendicular to 
the middle surface is ignored to improve numerical conditioning when the 
three dimensional element is employed. This constraint corresponds only 
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to a part of the usual assumptions of a two-dimensional shell theory. 
The relaxation of the requirement that straight lines perpendicular to 
the middle surface remain normal to the deformed middle surface permits 
the shell to experience shear deformation - an important feature in 
thick shell situations. 


Displacement Field in t he Element 

In the present study the current coordinates ‘x, are interpolated 

by the expression 


t x. = i + \ e 3i^ 

1 j-1 J 


and the displacement by 


V * £ + 7 ch j ^ e 3i ' e 31^ 

1 j-1 3 J 

n , i l . ,t+at:j _ tgj ii 

u. * e t + ? ch i ( e 3i e 3i' 

1 j-1 J J 


( 13 ) 


( 14 ) 

( 15 ) 


Here ‘u* and U J denote, respectively, the displacement and incremental 
displacement components in the x,-d 1 rection at the J-th node. The unit 
vectors \\ and e^ can be obtained from the relations 


• (f 2 x ‘e^/l^x 





( 16 ) 


where E 2 is the unit vector along the (global) x 2 -axis. If we assume 
that the angles e’ and s’ are very small, then we can write 
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Substituting Eq. (17) Into Eq. (15), we obtain 


u i ■ j, *j [u i + ? th j ( - + (i8a > 


or 

{u} = [T] {a} (18b) 

where {u} Is the column of three displacements at a point, {a} is the 
column of 5n (five per node) displacements: u^, 9 ^, e:j, j = 1,2,..,, n; i 
* 1,2,3, and [T] is the transformation matrix defined by Eq. (18a). 

Thus for each time step one can find the normal vectors from Eq. (16) 
and (17), and the incremental displacements at each point from Eq. (18) 
once the five generalized displacements at each node are known. 


Element Stiffness Matrix 

The strain-displacement equations (4) can be expressed in the 
operator form 


{e} = (A] {u q } (19) 

where {e} = {e n e 22 e 33 2e 12 2e 13 2e 23 J T , [A] is a function 
of t u 0 . j, and {u Q } is the vector of the components of the displacement 
gradient 


< u o> * t U l,l U 1 ,2 U 1 ,3 u 2,l u 2,2 u 2,3 u 3,l u 3,2 u 3,3^ (20) 

The vectors {u Q } and {e} are related to the displacement increments by 

{u Q } = [N] {u} = [N] (T] {a} (21) 

{e} = [A] [N] [T] {a} = [B]{a} (22) 

where [N] is the operator of differentials. 



Substitution of Eq. (22) into Eq. (8) yields 



where ‘(KJ, ‘lK NL l. (R}, and {F} are the linear and nonlinear stiffness 
matrices, force vector, and unbalanced force vectors: 


l lK, 1 = ; ‘lB! T lC1 t tBldV 
L V o 


^[K m1 1 - ; *IB! T IS] ^BldV 

NL w 

0 


(F) - J t IBl T (S(dV 0 
V o 


(24) 


Here [S] and {S} denote the matrix and vector, respectively, of the 
second Piola-Kirchhoff stress. 

Since we are dealing with laminated composite structures, the 
important thing is how to perform the integration through the 
thickness. One way is to pick Gaussian points through the thickness 
direction. This increases the computational time as the number of 
layers is increased, because the integration should be performed 
separately for each layer. An alternative way is to perform explicit 
integration through the thickness and reduce the problem to a two 
dimensional one. The Jacobian matrix. In general. Is a function 
of 5 1# t 2 , and c The terms in ; to the first power may be neglected, 
provided the thickness to curvature ratios are small. This 
approximation implies that derivative of x i with respect 
to 5 1# s 2 , and t are substantially the same at either end of a mid- 
surface-normal line. Thus the Jacobian (J1 becomes independent of ; and 
explicit integration can be employed. If c terms are retained in Ul. 
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Gaussian points through the thickness should be added. In the present 
study, it is assumed that the Jacobian is independent of c. 

Time Integration 

The Newmark Integration scheme Is used to convert the ordinary 
differential equations in time, Eq. (23), to algebraic equations. In 
the Newmark scheme, displacements and accelerations are approximated by 

t+4t {a} » *{ 4 ] + at 2 {i( + [(^ - B) t {a} + 6 t+4t {4(](»t) 2 

t+4 ‘{i} - ‘{i} + [(1 - ,)*{!} + T t+4t {i}l4t (25) 

where {a} is the generalized displacement vector of any point 
and 6 and y are the dimensionless parameters of the approximation. For 
the constant average acceleration case, we have s = ^ and y ~ and for 
the linear acceleration method e = g- and y = ^ (see [11]). 

Substituting Eq. (25) Into Eq. (23), and some algebraic 
manipulation leads to 

(a/lM] ♦ t (K]){a< k >} = t+at [R} - t+ 4 t {F< k - 1 >} + a 3 {P 4 } 

* a 2 |t(l> l ) • 4t (t(P 2! - t ( p 3 !)l (26) 

where 

a ° ' I^j 2 ’ = wt ’ a 3 * k - 1 ■ and 

[Ml - J 0. ‘[T] T ‘[T| dV 

\J ° ° 
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(Pl ] . J ‘W lT|dV 0 

0 

{P 2 }.J , 0 wt M (k *l>l T ,« 0 

V o 

{p,W, v »o ‘WITl*. 

0 

(?*} . J v *{i) [TldV 0 < 27 > 

0 

This completes the finite-element formulation of the 3-D degenerated 
element. 

nTSfUSSION OF THE NUM ERICAL RESULTS 
The results to be discussed are grouped into three major 
categories: (1) static bending, (2) natural vibration, and (3) 

transient response. All results, except for the vibrations, are 
presented in a graphical for*. All of the results presented here were 
obtained on an IBM 370/3081 computer with double precision arithmetic. 

Static Analysis 

Here we present a discussion of four example problems, all 
involving shell structures. 

1 Cylindrical Sh.il Subjected to Radia l Press_ure Consider a 
circular cylindrical panel of the type shown in Fig. 3. The shell is 
clamped along all four edges and subjected to uniform radial inward 
pressure. The loading is nonconservative, that is, the direction of the 
applied load is normal to the cylindrical surface at any time during the 
deformation. The geometric and material properties are 
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supported 
by a rigid 



diaphram: u = w= 0 



free edge 


Figure 3 Geometry of the cylindrical shell used in Problem 1 
of the static analysis. 



r . 2540 mm, a = b - 254 im, h * 3.175 mm, 

8 = 0.1 rad, E . 3.10275 kN/mm 2 , v = 0.3 

Due to the symmetry of the geometry and deformation, only one quarter of 

the panel Is analyzed. A load step of 0.5 KH/m 2 -as used in order to 

get a close representation of the deformation path. Fig. 4 contains the 

piot of central deflection versus the pressure. The solution agrees 

very closely with that obtained by Dhatt 113]. 

, n.thntrooic Cy »"dT Subjected tn internal Pressure Consider a 

clamped orthotropic (E 2 = 2 0 * 10 6 psi, * 3.75. G 12 /E 2 - 

0.625, v = 0.25) cylinder of radius R * 20" and length 20 
subjected to internal pressure, p 0 * 6.41/. psi. A mesh of 2x2 nine 
node elements is used to analyze the problem. The linear center 
deflections obtained by the 2-0 and 3-D elements are 0.0003764 in., and 
0.0003739 in., respectively. These values compare favorably with 
0.000366 in. of Rao 1141 and 0.000367 in. of Timoshenko's analytical 
solution 1151. The latter two solutions are based on the classical 

shell theory. 

In the large-deflection analysis the present results are compared 
with those of Reference 5. A value of 2.5 ksi Is used for the load 
step. Figure 5 contains a comparison of the present deflection wUh 
that of Reference 5. which used a 3-D degenerated element based on the 

updated Lagrangian approach. The agreement is very good. 

V nn._i,ver Cross- P ly fQ"/90"/0V. . . ) Spherical Shell Subject ed 

to uniform Loading Consider a spherical shell laminated of nine layers 

of graphite-epoxy material (E 1 /E 2 ■ 40, G i 2/ E 2 - °- 6, G 13 12 

G , V =.25 ), subjected to uniformly distributed loading, and simply 

i2 , . fr.ncifPrse deflection and tangential 

supported on all its edges (i.e., transverse oer .ect 
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Transverse load, 



l - 6 . 8 . 10 . 12 . 
Center deflection, w (in ran) 


Figure 4 


Load-deflection curve for the clamped cylindrical shell 




rotations are zero). A comparison of the load-def lection curves 
obtained by the present elements with those obtained by Noor [4] is 
presented- (for the parameters h/a = 0.01 and R/a = 10) in Fig. 6. The 
results agree very well with each other, the present 2-D results being 
closer to Noor's solution. This Is expected because Noor's element Is 
based on a shell theory. 

4. Two-Layer Cross-Ply and Angle-Ply (45 0 /-45 0 ) Shells Under 
Uniform Loading The geometry of the cylindrical shell used here is the 
same as that shown in Fig. 3. The shell is assumed to be simply 
supported on all edges. The material properties of individual lamina 
are the same as those used In Problem 3. A mesh of 2x2 nine-node 
elements in a quarter shell is used to model the problem. The results 
of the analysis are presented in the form of load-deflection curves in 
Fig. 7. From the results, one can conclude that the angle-ply shell is 
more stiffer than the cross-ply shell. 

The geometry and boundary conditions used for the spherical shells 
are the same as those used in Problem 3. The geometric parameters used 
are: R/a = 10, a/h = 100. The load-deflection curves for the cross-ply 

and angle-ply shells are shown in Fig. 8. From the plot it is apparent 
that, for the load range considered, the angle-ply shell, being stiffer, 
does not exhibit much geometric nonlinearity. The load-deflection curve 
of the cross-ply shell exhibits varying degree of nonlinearity with the 
load. For load values between 100 and 150, the shell becomes relatively 
more flexible. 

Natural Vibration of Cantilevered Twisted Plates 

Here we discuss the results obtained for natural frequencies of 
various twisted plates. This analysis was motivated by their relevance 
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Center deflection, w/h 






Nondi mens ionali zed center deflection, w/h 



Figure 7 Deflection versus the load parameter for two-layer 
composite cylindrical shell 



to natural vibrations of turbine blades. Consider 4n isotropic 
cylindrical panel with a twist angle e at the free end. Table 1 
contains the natural frequencies of a square plate for various values of 
the twist angle e and ratios of side to thickness. A 2x2 mesh and 4x4 
mesh of 9-node elements are employed to study the convergence trend. 

The results of the refined mesh are Included In the parentheses. The 
results obtained by using the 4x4 mesh are lower than those predicted by 
the 2x2 mesh, showing the convergence. The results agree with many 
others published in a recent NASA report. Table 2 contains natural 
frequencies of twisted plates for the aspect ratio of 3. 


Transient Analysis 

1. Spherical Cap Under Ax i symmetric Pressure Loading Consider a 
spherical cap, clamped on the boundary and subjected to axi symmetric 
pressure loading, p Q . The geometric and material properties are 
R = 22.27 in., h = 0.41 in., E = 10.5 x 10® psi, v = 0.3, 
o = 0.095 lb/in^, e = 26.67°, p Q = 100 psi. At = 10"® sec. 

This problem has been analyzed by Stricklin, et al. [16] using an 
axisynmietric shell element. In the present study the spherical cap is 
discretized into five nine-node 2-D and 3-D elements. Figure 9 contains 
the plot of center deflection versus time. The present solutions 


obtained using the 3-D and 2-D elements are in excellent agreement in 
most places with that of Stricklin et al [16]. The difference between 
the solutions is mostly in the regions of local minimum and maximum. 

2. Two-Layer Cross-Ply Plate Under Uniform Load A cylindrical 
shell with a = b = 5", R = 10", h = 0.1" is simply-supported on the four 
edges, is analyzed. The shell is laminated by 2 layers (0°/90°) and 

- A 

exerted by a uniform step load P = = 50. Figure 10 contains a plot 

E 2 »T 
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* 2x2, 9-node mesh 
**3x3, 9-node mesh 
+ 4x4, 9-node mesh 
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Table 2 Natural Frequencies of Twisted Rectangular Plates 
(b/a = 3, 3x3 mesh of nine-node elements) 

~ = ajb 2 /pT?U , D = — Ll? — , v = 0.3 

12(l-v) Z 


a 

h 

Twi st 
Angl e 

1 

2 

3 

Mode 

4 

5 

6 

7 


0° 

3.4150 

20.8772 

21 .6190 

65.9706 

66.2590 


1 27 .256 


15° 

3.4009 

20.8798 

22.1 118 

21.5032 

68.0938 

69.3258 

130.284 

20 

30° 

3.3598 

19.4048 

25.3743 

60.2183 

73.5180 

77.4493 

138.176 


45° 

3.2956 

17.5289 

29.8404 

58.2600 

80.9488 

88.5245 

148.8975 


60° 

3.2136 

15.7431 

34.8827 

55.8921 

89.2028 

100.7760 

155.070 


0° 

3.3908 

15. 

551 

19.124 

21.065 

59.924 

61.949 


15° 

3.3161 

15. 

192 

19.231 

21.572 

60.088 

60.830 

5 

30° 

3.3336 

14. 

379 

19.549 

22.81 1 

60.576 

58.472 


45° 

3.2674 

13. 

449 

20.060 

24.404 

51 .360 

55.874 


60° 

3.1833 

12. 

548 

20.741 

26.139 

62.416 

53.381 
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Center displacement, wxlO (in inch) 



-1 .0 J 


Figure 9 Center transverse displacement versus time 

for a spherical cap under axi symmetric dynamic 
loading ( load * 100 psi.) 
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Center deflection, wxlO 



Figure 10 Center deflection versus time fo 
cross-ply cylindrical shell subj 
uniform step load 


Of the center defection versus time for 2-0 end 3-0 elements. The time 
step used is at - 0.1 x 10‘ 4 sec. The solutions obt.ined using the two 

elements are in good agreement. 

3 T ■ Anglej^ f«V-45-» Spherical then Under Uniform 

Loading Consider a spherical shell with » ■ b • 'O’. R ' 20 “ and " ’ 

0.1" f simply supported at four edges and Is everted by a uniform step 

load. The shell consists of two layers. («•/-«•>. Figure 11 contains 

the plot of center deflection versus time for P -50 and P - 500 with 

time step 0.2 x W* sec. For the small load the curve is relatively 

smooth compared to that of the larger load. This is due to the fact 

p - 50 is smaller compared 

that the geometric nonlinearity exhibited a 
to that at P * 500. CONCLUSIONS 

The present 3-0 degenerated element has computational simplicity 
over a fully three-dimensional element, such as those developed in 1171. 
and the element accounts for full geometric nonlinearities in contrast 
t „ 2-0 elements based on shell theories. As demonstrated via numerical 
examples, the deflections obtained by the 2-0 shell element deviate from 
those obtained by the 3-D element for deep shells. Further, the 3-0 
element can be used to model general shells that are not necessary 
doubly-curved. For example, the vibration of twisted plates cannot be 
studied using the 2-0 shell element discussed in [121. Of course, the 
3-0 degenerated element Is computationally more demanding than the 2 0 
shell theory element for a given problem. In su«ary. the present 3-0 
element is an efficient element for the analysis of laminated composite 
plates and shells undergoing large displacements and transient motion. 

The 3-0 element presented herein can be modified to include thermal 
stress analysis capability and material nonlinearities. While the 
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Center deflection 




inclusion of thermal stresses is a simple exercise, the inclusion of 
nonlinear material effects is a difficult task (see (18-20]). An 
acceptable material model should be a generalization of Ramberg-Osgood 
relation to an anisotropic medium. Another area that requires further 
study is the inclusion of damping effects, which are more significant 
than the shear deformation effects. 
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PART 3 


NONLINEAR MATERIAL MODELS FOR COMPOSITE PLATES AND SHELLS 

K. Chandrashekhara and J. N. Reddy 
Department of Engineering Science and Mechanics 


SUMMARY 

Nonlinear material models for laminated structures are described 
and their Incorporation In the finite-element formulation of laminated 
plates and shells Is presented. Numerical results for several sample 
problems of plates and shells are presented and validated by comparison 
with those available In the literature. 

INTRODUCTION 

Composite materials are known to exhibit significant non- 
linearities in stress-strain behaviour even at low strains. Most of the 
currently used matrix materials In composites have high strain 
capabilities and the Investigation of the bending of composite shells 
undergoing large deformation, yielding is apt to occur and Its effect 
must be accounted for In the analysis. The nonlinearity Is not 
Isotropic but varies with direction, as do the elastic properties. 

Models for such elastic-plastic behavior of orthotropic and anisotropic 
materials are not well developed. 

The total stress-strain laws are mathematically more convenient 
than incremental laws but are physically not sound. The criterion 
approximately describing the yielding of Isotropic material Is that of 
von-MIses. The simplest yield criterion for anisotropic material Is 
therefore one which reduces to von-Mises law when the anisotropy is 
vanishingly small. Hill's yield criteria assumes relatively simple ease 
of orthotropic anisotropy, that Is, there are three mutually orthogonal 
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planes of symmetry at every point and the Intersection of these planes 
are considered as the principal axes of anisotropy. Fiber reinforced 
composite structures almost Invariably possess this kind of symmetry. 

In the present study a nonlinear material model Is developed for 
composite plates and shells, and numerical results for bending are 
presented using the finite element method as exact solutions are not 
tractable for elastic-plastic problems Involving complex geometries. 


MATERIAL MODEL 

In the present model. Hill's anisotropic yield criteria for 
elsatic-perfectly plastic material is used. Hill's [1] yield function 


f(o 1;J ) - F(o 2 - o 3 ) 2 + G(a 3 - a x ) 2 + H( 0l - o 2 ) 2 

+ 2Lo 2 3 + 2Mo 2 3 + 2No 2 2 = 1 (1) 

where F, G, H, L, M, N are parameters characteristic of the current 


state of anisotropy given by, 



and X, Y, Z are the tensile yield stresses In the principal direction of 
anisotropy and R. S, T. are the yield stresses In shear with respect to 

the principal axes of anisotropy. 

It should be noted that Hill's criteria is based on the assumption 
that the superposition of a hydrostatic stress does not Influence 
yielding and there is no Bauschinger effect. Also, the yield criterion 
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has this form when the principal axes of anisotropy are the axes of 
references. 

For a plane stress state in the 1-2 plane with transverse shear, 
equation (1) reduces to: 

f * (G + H)o 2 + (F + H)o 2 - 2Ha l a 2 

+ 2lo\ z + 2Mcj 3 + 2Hc\ z » 1 (2) 

For an isotropic material; 


V 

the yield stress In uniaxial tension and according to the von-Mises 
yield criteria [2] 


R * S = T = -S 
✓3 

Therefore, F = G = H = — ^ and 21 = 2M = 2N = -^ and equation (2) 

2c* a o 

becomes , ° ° 

- 2,2 2 2 2 v 2 

f “ °1 °2 " °1°2 3 * ^°23 °13 °12^ _ °o 

which Is the familiar von-Mises yield criteria. 

If the principal axes of anisotropy 1,2 do not coincide withthe 

reference axes x, y, but are rotated by an angle 6, then the stresses in 

equation (2) are obtained using the transformation as: 

2 2 

o. = o cos e + o sin e + o sine cose 
ix y xy 

2 2 

o ? = o sin e + a cos e - a sine cose 

c x y xy 


23 


-°X2 Sl " 9 


o yz cose 


o ^3 * o x2 cose + o y2 sine 


12 


= -2a sine cose + 2a sine cose + a (cos^e - sin^e) 
x y xy 


Elastic-Plastic Constitutive Equations 

In the Incremental theory of plasticity, the total strain increment 
is the sum of the elastic and plastic components 


76 



de 3 dc e + de p 


(3) 


The elastic strain increment is related to the stress increment by 
Hooke's law as, 

de e * [D e ] -1 do 

where [D e ] is the elastic modulus matrix which for orthotropic material 
takes the form. 


[0 e l = 


1 " v 12 v 21 

v 12 E 2 

l' v 12 v 21 


0 

0 

0 


v 12 £ 2 

1-v 12 v 21 


1 " v 12 v 21 

0 

0 

0 


0 

0 

G 23 

0 

0 


0 

0 


*13 


0 

0 

0 
0 

G 12j 


(5) 


The normality rule for an associated plastic flow is. 


<*' p - <& & 

where dx Is the positive proportionality constant, evaluated using the 
condition that during the plastic deformation, the stresses remain on 
the yield surface so that, 

df 3 — da 3 0 
3o 

The stress-strain relation in the plastic range is given by [3], 

do 3 [D ep ]dc 


where 


l° ep l ■ to”] - TifTT 


t° e H£H£)Vl 


i£n° , 0 


( 6 ) 
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Hence the modification called for In the elastic-plastic analysis 
would be solely the replacement of the elasticity matrix [D e ] by the 
elastic-plastic matrix [D ep ] for the yielded elements at the successive 
stages of calculation. It should be noted that the [D ep ] matrix Is 
populated and accordingly the transformation of the stress-strain 
relation from the material axes, ({o} * [0 ep ]j 2 {e}), to the shell 
coordinate axes, ({0} = [D ep ] {c} ) , will be modified as shown In 
Appendix I. 

FINITE ELEMENT FORMULATION 

Consider a laminated shell constructed of a finite number of 
uniform thickness orthotropic layers, oriented arbitrarily with respect 
to the shell coordinates The orthogonal curvilinear 

coordinate system (^,^.0 is chosen such that and $ 2 - curves are 
lines of curvature on the midsurface c=0, and c-curves are straight 
lines perpendicular to the surface c=0. 

For the small displacement Sanders shell theory which accounts for 
transverse shear deformation, the strain displacement relations are 
given by [ 5] , 

e 1 = e i + C *1 

where 

e l 3X^ + * 1 3Xj 

0 + . _!!i 

c 2 * 3Xj_ R 2 * k 2 3X 2 

0 . !!i + mi . . !!i + 

e 6 ~ 3X 2 ’ *6 " 3X 2 3Xj C 0 '3X^ 3X 2 ^ 

0 A + ^ ^ 

e 4 " *2 3 x 2 ' R 2 
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o + IHi Hi 

e 5 “ *1 3X X ' R x 

c = — (i — 4-) » dx. * «4<U, ^ 

c o 2 'Rj R 2 1 11 

Here R 1 (1 * 1.2) are the principal radii of curvature, u 1 are the 
displacements of the reference surface along 5,(53 - t) ««, *, and * 2 
are the rotations of the transverse normals about the 5 2 and 5,-axes 

respectively. 

The stress-strain relations, transformed to the shell coordinates, 
are of the form 

{0} * [ Q 1 { e } 

where are the material properties of k th -layer (see Appendix I). 
The principle of virtual work for the present problem is given by 


0 = Z [J + <4 k) Se 2 + 6eg + oJ k) 6 « 

1=1 0 


+ oi k ^ «cc - qau^ajOgd^d^ldt 


(7a) 


= J [HjdeJ + N 2 6e° + N 6 5 c° + M 1 fi< 1 + + * 6 &< 6 

a 

+ N 4 5e° + N 5 6e° - q^Ja^d^d^ (7b) 

where q is the distributed transverse load, H i and M i are the stresses 
and moment resultants. 

L c k 

(N . ,M.) = l J 0^(1, c)dc (1 = 1,2, 6, 4, 5) 

1 1 k=1 H-i 

Here c k ) are the c-coordinates of the k th layer, and L Is the 

total number of layers in the laminated shell. 
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It should be noted that the equations of equilibrium can be derived 
from Eq. (7b) by Integrating the displacement gradients In by parts 

and setting the coefficients of «u i to zero separately. We obtain 

3N 1 3 N 5 

1 + jfc (N « + C A) - B 5 ' 0 


ax 


1 


'o' ’6' 
aN, 


a 

axi 

< N 6 - 

e oV + 

c 

3X 2 

+ TT- 0 
K 2 

aN 5 

ax^ 

+ 2!k 

8 X 2 

M 1 

** 
*2 ' 

. q) - 0 


aM 1 

aM, 

+ — - - 

N c * 

0 


axj 

a x ^ 

5 



aMg 

aM 2 

+ ” 

H. * 

0 


axj 


4 



The resultants (N^, M^) are related to (e°,ic..) by, 
o 

t- D,_<. 

1.J - T.2.6,4,5 

a ,p = 1,2,6 (with 1*1 for i = 1,2,6) 


N 


1 


A 1j e j + B -< 


ip p 


M = 8 .c, + 0* 
a aj j ip p 


( 8 ) 


Here A^j, and D^- denote the extenslonal, flexural-extensional 
coupling, and flexural stiffnesses of the laminate: 


1 J (l,c,c 2 )di: 

"- 1 H-i J 


( A i j ,B i j ,B i j) 

In the unabridged notation equation (8) takes the form: 


(9) 
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v 


_A 11 a 12 A 16 -14 -15 B 11 B 12 B 16 


0 1 
e l 


n 2 


a 12 a 22 a 26 a 24 a 25 b 12 b 22 b 26 


u 

e 2 

A 


N 6 


a 16 A 26 A 66 -46 -56 B 16 B 26 B 66 


U 

e 6 

A 


N 4 


a 14 A 24 -46 A 44 A 45 Sl4 ^24 &46 


U 

e 4 

_ f 

< 

N 5 

*x 

a 15 A 25 -56 A 45 A 55 §15 ^25 §56 

n 

«! 


Mi 


B 11 B 12 B 16 -14 -15 °11 °12 °16 


*1 


M 2 


B 12 B 22 B 26 B 24 B 25 D n D 22 0 26 


< z 


• M 6. 


b 16 b 26 b 66 b 46 b 56 d 16 o 26 d 66 __ 


k 6 

The underscored coefficients are due to material nonlinear stri 


( 10 ) 


relationship. It should be noted that the coefficients A 44 , A 45 and A 55 
defined in equation (9) has to be corrected for the parabolic variation 

of the transverse shear stress, as 


( A 44’ A 45 ,A 55^ 


L 

i I 

k=l c 


c k 

k-1 


fk 2 0 (k) 

(*r44 • 


k k 0^ k 2 0^ )dc 
k l k 2 Q 45 * 2 g 55 } ; 


( 11 ) 


where k are the shear correction factor. 

A typical finite element Is a doubly-curved shell element whose 
projection is an isoparametric rectangular element. Over the typical 
shell ii (e) . the displacements Interpolated by 

expressions of the form, 


ii j 

u. = E Uai>e(x. ,x ? ) , i - 1*2,3 

1 j-1 1 J 1 


*1 j 

= E *j4».(x.,x ? ) , i * 1,2 

i j =1 1 J L 


( 12 ) 


where ip 


are the interpolation functions, 

j 


and u^ and 4>'- 


are the nodal 
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values of u^ and $.j, respectively. For a nine node quadratic element 
the element stiffness matrix is of order 45x45. 

Substitution of equation (12) into the virtual work principle, Eq. 
(7b) yields an element equation of the form 

(KIM - {F} (13) 

where {a} * {{uj}, {u 2 }, {u 3 }, {♦j}, {* 2 }} T , (KJ Is the element 
stiffness and {F} Is the force vector. In the Interest of brevity, the 
coefficients of stiffness matrices are Included In Appendix II. 

It should be noted that the underscored coefficients in Eq. (10) 
are also redefined like the shear coefficients in Eq. (11) and reduced 
integration Is performed for the terms arising in the element stiffness 
matrices due to the presence of these coefficients to avoid the so- 
called locking effect. 


NUMERICAL RESULTS 
The Parameters of Anisotropy 

When considering the modeling of a material system, one must always 
survey the availability of material property data. In the present 
theory, to describe fully the state of anisotropy, the six independent 
yield stresses in Hill's criteria are needed to be known from uniaxial 
tests. For numerical results, two typical composite materials namely, 
boron/epoxy and graphite/epoxy are considered with the following 
material constants: 

Boron/Epoxy 

E x * 30.0 x 10 6 psi , E 2 = 3.2 x 10 6 psi 
Gj 2 = 1.05 x 10 6 psi , v^ 2 * 0.21 , G 2 2 * Gj 3 = G 12 
X = 195 x 10 3 psi ; Y = Z * 12.5 x 10 3 psi 
R = S = T = 18.0 x 10 3 psi 
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Graphite /Epoxy 

E^ = 18.88 x 10® psi ; E 2 * 1.376 x 10® psi 
G 12 « 0.688 x 10 6 psi ; v 12 =» 0.343 ; G 23 * G 13 

X * 222.7 x 10 3 psi ; Y = Z » 6.35 x 10 3 psi 
R = S = T = 9.92 x 10 3 psi 


G 


12 


Solution Procedure 

The solution of the elastic plastic problem Is reached by an 
Incremental and Iterative procedure. The direct Iteration technique Is 
followed in the present analysis. 

For each load increment, the system of equations are established by 
assembling the element matrices and the displacement {a} is obtained 
from Eq. (13) . Consequently, the state of stress and the value 
of f(o. .) are calculated for each element. If f < 0, then the process 
is elastic and the material matrix is obtained from equation (5). If f 
> 0, then the total stresses are readjusted so as to make f « 0 and the 
elastic-plastic matrix is calculated from Eq. (6). Once the convergence 
is achieved, the next load increment is applied and the iteration 
procedure is repeated. 

If the application of a small load increment causes very large 
deflection, the calculation is stopped and the limit load is considered 
to be found. 

Sample Problems 

The present elastic-perfectly plastic formulation is applied to a 
variety of bending problems using 2x2 mesh of a nine noded quadratic 
element. The shear correction factors k* = k* were taken to be 5/6. 
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All computations were made using an IBM 3081 processor with double 
precision arithmetic. 

The results of the sample problems are presented and compared, if 
possible, with the existing solutions to evaluate the present 
formulation. 

1. Cylindrical Shell Roof A cylindrical shell subjected to 
uniform vertical loading Is considered. Due to symmetry, only a 
quadrant of the shell was analyzed. The geometry and modeling of the 
shell roof are shown in Fig. 1. The material behaviour is studied with 
the properties: 

E : = E 2 = 2.1 x 10 4 MN/m 2 ; v = 0.0; 

G 12 = 1.05 x 10 4 MN/m 2 ; G 23 = G 13 = G 12 

X = Y = Z = 4.1 Mn/m 2 ; R = S = T = 2.367 MN/m 2 

The results obtained for the vertical displacement at the central 
point of the free edge A versus loading was shown in Fig. 1. The 
solution obtained compares well with those reported In Ref. [6]. The 
apparent discrepancy can possibly be due to a different boundary 
condition on the curved edges and the type of material model used. 

2. Simply-Supported Square Plate A uniformly loaded simply 
supported square plate was studied In the second example. The geometry 
of the plate Is shown in Fig. 2. The following material properties were 
considered: 

E. = E 2 * 10 x 10 6 psl ; v = 0.3 
G 12 = 3.846 x 10 6 psl ; G 23 * G 13 = G 12 
X = Y = Z = 144,000 psl ; R * S = T = 83,138.4 psi 
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Figure 1. Load-deflection curves for a cylindrical panel 
under uniform transverse load 





A non-dlmenslonalized plot of the centre displacement of the plate 
versus the load are shown In Fig. 2. The results are compared with 
those presented in Ref. 17]. 

V Two Laver Cross-Ply [0/90] and Angle-Ply f-45/45] Simply 
Supported Spherical Shells Figure 3 contains the results for the cross- 
ply shell made of two typical materials, namely, boron/epoxy and 
graphite/epoxy under uniform load. For a given load, the shell made of 
graphite/epoxy deflects more than the shell made of boron/epoxy which Is 
stlffer, but experiences small degree of nonlinearity. 

Figure 4 shows nonlinearity exhibited by the graphite/epoxy cross- 
ply and angle-ply shells under uniform load. Clearly, the angle-ply 
shows greater displacement and also nonlinearity than the cross-ply for 

the same load. 

Figure 5 shows the material behaviour for the boron/epoxy cross-ply 
shell under concentrated load. 

a. Clamped Cylindrical Cross-Ply ( 0/901 Shell Under Uniform Load 
The geometry of the shell is shown in Fig. 6. The shell is made of 
grpahite/epoxy and the plot of displacement versus load are shown in 

Fig. 6. 

CONCLUSIONS 

A finite element model based on Sander's shell theory, accounting 
for the transverse shear strains is used for the elastic-plastic 
analysis of lamianted composite shells. The parameters of anisotropy 
reflect the plastic material response by correcting the stress 
components in the Hill's yield function. Numerical results are 
presented for isotropic and laminated shell of cylindrical and spherical 
geometry to demonstrate the validity and efficiency of the present 
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Figure 3. Load-deflection curves for a simply supported 
spherical shell under uniform transverse load 
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\k 





Fiqure 5. Load-deflection curves for a simply supported spherical shell 
under point load at the center (see Fiqure 3 for the qeonetry) 



100*1 Linear- 
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approach. For the Isotropic case, the present results are in good 
agreement with those available in the literature. 
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APPENDIX I 


Transformation of the stress-str ain m atrix In Equation .£61 
Let the elastic matrix In the material axes (1,2) be (D p li2 and 
the body axes be [O ep l x y 


tD ep J lz - IC1 


[D ep ] xy = [Q] 


Cll c 12 C 13 C 14 C 15 

°22 C 24 C 25 

c 33 c 34 ^35 

\ 

sym C45 

c 55_ 

Qll Q 12 Q 13 hi 

S Q 22 Q 23 Q 24 ^25 

Q 33 Q 34 ^35 

S 

sym Q44 Q45 

q 55_ 


then the transformation 14] Is given as, (with m = cose, n • sine) 

Q n - m 4 C u + 2m 2 n 2 (C 12 * 2C 33 ) - 4mn(m 2 C n * n 2 C 23 ) ♦ 

«12 * » 2 " 2 (C U * C 22 - 4C 33 ) * 2mn(m 2 - n 2 )(Cj 3 - C 23 ) * 0" 4 * " 4 > C 12 
q 13 . m 2 (m 2 - 3n 2 )C 13 ♦ mn(m 2 C 11 - n 2 C ?2 - 0» 2 - " 2 )(C 12 * 2C 66 ) 1 

+ n 2 (m 2 - n2 )C 2 g 

Q 14 = m 3 C u - mn((2C^ - C 15 )m - (C ?4 - 2C 35 )n) + C^n 3 
Q 15 = m 3 C 15 - mn[(C 14 + 2C 35 )m - (C 25 + 2C 34 )n] - C 24 n 
q 22 = n 4 C u + 2m 2 n 2 (C 12 + 2C 33 ) + 4mn(m 2 C-, 3 ♦ n 2 C 13 ) * m 4 C 22 
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Q 23 = m 2 (m 2 - 3n 2 )C 23 + mn[n 2 C u - m 2 C 22 + (m 2 - n 2 )(C 12 + 2C 33 ) ] 

+ n 2 (3m 2 - n 2 )C^ 3 

Q 2 4 * >" 3 C 24 + mn[(C 25 + 2C 34 )m + (C 14 + 2C 35 )n] + C 15 n 3 

Q 2 5 * m C 2 g - mn[(C 2 ^ - 2C 3 g)m - (C^g - 2C 3 ^)n] - Cj^n^ 

Q33 * m n (C^ + C 22 - 2 Cj 2 ) - 2mn(m 2 - n 2 )(C 22 - C 33 ) + (m 2 - n 2 ) 2 C 33 

Q34 * ^ C 34 + nC 35 )(m2 ~ p2 > * m2n ( C i4 - c 24 ) + nrn 2 (C ls - C 25 ) 

Q 35 * (mC 35 - nC 34 )(m 2 - n 2 ) ♦ m 2 n(C ls - Cgg) ♦ mn 2 (C 24 - C 14 ) 

^44 = m C 44 + 2wnC 45 + n2c 55 

Q45 “ n ^45 " mn (^ 44 ~ ^55) 

^55 = m2c 55 " 2mnC 45 + p2c 44 

The underscored terms are due to material nonlinearity for an 
orthotropic material. Also note that the constitutive matrix is no 
longer orthotropic. 
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APPENDIX II 


Stiffness Coefficients 




|K 14 | . B n lS U l * B 16 ([S 12 ] ♦ IS 21 )) * B 66 IS 22 1 * 
C o (0 16 ls21 ' + °66 Is22|) " *55 IS001 

- A 15 tS 10 l * A 56 IS ; °] - fr (B 15 iS 01 i » B 56 IS° 2 n 

♦ W* 20 ' 




[K 15 ] = B 12 [S 12 ] + B 16 [S U ] + B 26 (S 22 1 + B 66 [S 21 ] 

* C o(0 26 [S 22 I + D 66 IS 21 ]) -^ a 45 is 00 i 

* A 1 4 |s1 ° 1 + A 46 ls20) - EJ < B 25 |S ° 21 + B 56 IS ° l » 

* C 0 B 46 [S2 °I 

[K 21 ] = [K 12 ] T 

[K 22 ] = A 22 [S 22 ] + A 26 ([S 12 ] + [S 21 ]) + A 66 [S U ] - 2C 0 B 66 [S U ] 

- c o (b 26 ([s 12 ] + [s 21 ]) - c 0 d 66 [s u ]) - ^ [S 00 ] 

- jq ( A 24< ts2Ql + lS ° 2) ) + A 46 (lsl0] + [s01] + C 0 IS 10 ]) 

- 

IK 23 ] - jq- (AjjIS 20 ] + A 16 1S 10 ]) + i- (A 22 [S 20 | + A 26 [S 10 1) 

- c o<sf + V' [S1 °' * k ( V S ° 2 1 + V s01 D 

+ A 24 ([S 22 ] - ij I S 00 1 ) + a 25 is 21 i + a 46 is 12 | 
^2 

* V s11 ' - e^e; t s °°) - W s ' 2 ' - W* 11 ' 

IK 24 ] = B 12 [S 21 ] + B 26 [S 22 ] + B 16 (S U ] + B 66 [S 12 ] 

- C 0<°16' SU 1 * °66 IS12 D - EJ A 45f S °°l 

+ A 25 |S201 * A 56 |S101 * r 2 < B 14* s01 ' + 8 46 [S ° 2, > 
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[n 25 l . b 22 Is 22 ] + B 26 (|S 21 1 + is 12 ]) ♦ B 66 IS U ] 

- C o (D 26 lslZ| * D 66 lsll|) ' *1 A 44 tS °°' 

♦ A 24 [S 2 °| ♦ A a6 is 10 i - (B 24 (S 02 1 ♦ b 46 1S 01 ]) 

- wfl!i 

[K 31 ] « [K 13 ] T 
IK 32 ] - [K 23 1 T 
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,35 


IK"1 * A._tS 10 ] + A..[S"1 


. 20 , 


45 


44 1 


+ (V + V )[s01l + ( T + ? )|s ° 21 

+ (lf + Rf)l S ° 0 l + B 24 | S 22 J + B. i: rs 21 | 


' 46 ' 



* B 2S [S 12 1 

* a S6 fsll l 



[K 41 ] 

- [K 14 ] T 




[1C 42 ) 

• [k 24 i t 




[K 43 ] 

= [K 34 ] T 




IK 44 ] 

■ 0 U [S U 1 ♦ d 16 ([s 12 j 

+ IS 21 ]) 

+ °66 |S2? 1 + 4 55 IS °°I 


- b 15 (is 10 : 

1 + [$ 01 ]) + 

B 56 ((S 20 

1 + [S 02 ]) 

[K 45 ] 

■ 0 12 (S 12 ] ♦ 0 16 [S U ] + 

d 26 [s 22 i 

♦ B 6 6 [s21 ' + A 45 [S °° 


+ B 14 [S 10 | 

* B 46 ts 20 i 

+ B 25' S ° 2 

1 + b 56 is 01 i 

[K 51 ] 

- ( K 1S ] T 




[K 52 ] 

- [K 25 ) T 




[K 53 ] 

- [K 35 ] 1 




[K 54 ] 

■ (k 45 ! t 




[K 55 ] = 

d 22 [s 22 ] + 

d 26 ([S 12 ] * 

IS 21 ]) + 

D 66 1s11 ' + A 44 1 B°° 1 


+ B 24 ((S 20 | 

+ [S 02 ]) + 

B 46([S 10 i 

1 + IS 01 ]) 


A 
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where 


and the 


3^.: 00 r 

joS s r — 1 — 1 dx.dx,, , S . . = j 

Mi J ax 3X. 12 'J 


'ij " J e 3x a 3x s 

Q 


underscored terms are due to material 


*1<»jd x id x 2 

nonlinearity. 
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